Streamer Branching rationalized by Conformal Mapping Techniques 
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Spontaneous branching of discharge channels is frequently observed, but not well understood. We 
recently proposed a new branching mechanism based on simulations of a simple continuous discharge 
model in high fields. We here present analytical results for such streamers in the Lozansky-Firsov 
limit where they can be modelled as moving equipotential ionization fronts. This model can be 
analyzed by conformal mapping techniques which allow the reduction of the dynamical problem to 
finite sets of nonlinear ordinary differential equations. Our solutions illustrate that branching is 
generic for the intricate head dynamics of streamers in the Lozansky-Firsov-limit. 
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When non-ionized matter is suddenly exposed to 
strong fields, ionized regions can grow in the form of 
streamers. These are ionized and electrically screened 
channels with rapidly propagating tips. The tip region is 
a very active impact ionization region due to the self- 
generated local field enhancement. Streamers appear 
in early stages of atmospheric discharges like sparks or 
sprite discharges [1,2], they also play a prominent role in 
numerous technical processes. It is commonly observed 
that streamers branch spontaneously [3,4]. But how this 
branching is precisely determined by the underlying dis- 
charge physics, is essentially not known. In recent work 
[5,6], we have suggested a branching mechanism from first 
principles. This work drew some attention [7,8], since 
the proposed mechanism yields quantitative predictions 
for specific parameters, and since it is qualitatively dif- 
ferent from the older branching concept of the "dielec- 
tric breakdown model" [9-11]. This older concept actu- 
ally can be traced back to concepts of rare long-ranged 
(and hence stochastic) photo-ionization events probably 
first suggested in 1939 by Raether [12]. Therefore, it 
came as a surprise that we predicted streamer branch- 
ing in a fully deterministic model. Since our evidence for 
the phenomenon was mainly from numerical solutions to- 
gether with a physical interpretation, the accuracy of our 
numerical scheme was challenged [13,14]. Furthermore, 
some authors have argued previously [15,16] that in a de- 
terministic discharge model like ours, an initially convex 
streamer head never could become locally concave, and 
that hence the consecutive branching of the discharge 
channel would be unphysical. 

Therefore in the present paper, we investigate the is- 
sue by analytical means. We show that the convex-to- 
concave evolution of the streamer head with successive 
branching is generic for streamers in the Lozansky-Firsov 
limit [5,6,17]. We define the Lozansky-Firsov-limit as 
the stage of evolution where the streamer head is al- 
most equipotential and surrounded by a thin electro- 
static screening layer. While in the original article [17], 
only simple steady state solutions with parabolic head 
shape are discussed, we will show here that a streamer 
in the Lozansky-Firsov limit actually can exhibit a very 



rich head dynamics that includes spontaneous branch- 
ing. Furthermore, our analytical solutions disprove the 
reasoning of [15] by explicit counterexamples. Our ana- 
lytical methods are adapted from two fluid flow in Hclc- 
Shaw cells [18-22]. But our explicit and exact solutions 
that amount to the evolution of "bubbles" in a dipolc 
field [23], have not been reported in the hydrodynamics 
literature either. 

The relation between our previous numerical investi- 
gations [5,6] and our present analytical model is laid in 
two steps. First, numerical solutions show essentially the 
same evolution in the purely two-dimensional case as in 
the three-dimensional case with assumed cylinder geome- 
try [5,6]. Because there is an elegant analytical approach, 
we focus on the two-dimensional case. This has the ad- 
ditional advantage that such two-dimensional solutions 
rather directly apply to, e.g., discharges in Corbino discs 
[24]. Second, we use the following simplifying approxi- 
mations for a Lozansky-Firsov streamer: (1) the interior 
of the streamer is electrically completely screened; hence 
the electric potential ip is constant; hence the ionization 
front coincides with an equipotential line, (2) the width 
of the screening layer around the ionized body is much 
smaller than all other relevant length scales and in the 
present study it is actually neglected, (3) the velocity of 
the ionization front v is determined by the local elec- 
tric field; in the simplest case to be investigated here, 
it is simply taken to be proportional to the field at the 
boundary v = c \7ip with some constant c (for the va- 
lidity of the approximation, cf. [25,26]). Together with 
~S/ 2 ip — in the non- ionized outer region and with fixed 
limiting values of the potential tp far from the streamer, 
this defines a moving boundary problem for the interface 
between ionized and non-ionized region. We assume the 
field far from the streamer to be constant as in our sim- 
ulations [6]. Such a constant far field can be mimicked 
by placing the streamer between the two poles of an elec- 
tric dipole where the distance between the poles is much 
larger than the size of the streamer. 

When the electric field points into the x direction 
and y parametrizes the transversal direction, our two- 
dimensional Lozansky-Firsov streamer in free flight in a 
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homogeneous electric field is approximated by: 



vV(z,y) = o, 

-Vtp(x,y) -> E x, 
tp(x,y) = 0, 

Vbound = C VwL , , 
~ Ibound. ' 



outside the streamer, (1) 

far outside the streamer, (2) 

inside the streamer, (3) 

velocity of the boundary, (4) 



where x is the unit vector in the x direction, and we 
have chosen the gauge such that the potential inside the 
streamer vanishes. The asymptote (2) implies that the 
total charge on the streamer vanishes; otherwise a con- 
tribution oc 1/ ' \Jx 2 + y 2 has to be added on the r.h.s. of 
Eq. (2). 

Similar moving boundary problems arise in Hele-Shaw 
flow of two fluids with a large viscosity contrast [18,19]: 
Lozansky-Firsov streamers and viscous fingers on the 
present level of description can be identified immediately 
by equating the electric potential ip with the fluid pres- 
sure p [25,26]. To such problems, powerful conformal 
mapping methods [20-22] can be applied. Most work 
with this method is concerned with viscous fingers in a 
channel geometry, i.e., with boundary conditions on a 
lateral external boundary that cannot be realized in an 
electric system. A few authors also study air bubbles 
within a viscous fluid, or viscous droplets in air, mostly 
under the action of flow fields generated by one source 
or one sink of pressure, i.e., by monopoles. On the other 
hand, the approximation (l)-(4) describes streamers in 
free space between two electrodes as in [6]. With the 
asymptote (2), this is mathematically equivalent to air 
bubbles in a dipole field. This case has not been studied 
in detail. It is known that any ellipse with the main axes 
oriented parallel and perpendicular to the direction of the 
dipole is a uniformly translating solution of this problem 
[27]. The time dependent solutions of [28] do not apply 
to streamers since the boundary condition on the mov- 
ing interface is different. Rcfs. [23] and [29] study how 
and when cusps in the interfaces of droplets and bubbles 
emerge when these are driven by multipole fields. But 
for bubbles in a dipole field, again only the steady state 
ellipse solutions are given [23]. 

In the present paper, we therefore apply conformal 
mapping methods to the evolution of "bubbles" in a 
dipole field in a Hele-Shaw experiment and proceed be- 
yond the steady state ellipse solutions. We identify the 
general structure of time dependent solutions of (1)- 
(4) . The analytically derived exact solutions show how a 
streamer head can become flatter, concave and branch as 
observed numerically [5,6]. Rather than a pole decompo- 
sition [20] , we derive a decomposition into Fourier modes 
of the circle and calculate an equation for the non-linear 
dynamical coupling of their amplitudes. 

In detail, this is done in the following steps: 
(i) The spatial coordinates are expressed by the complex 
coordinate z = x + iy. According to standard complex 
analysis, finding a real harmonic function <p(x, y) solving 



the Laplace equation (1) in a given region is equivalent 
to finding a complex function $>(z) that is analytical in 
the same region and has real part Re <&(z) = tp(x, y). 
(ii) A conformal map from the interior of the unit circle 
to the exterior of the streamer or "bubble" is constructed. 
Including the point at infinity, the region outside the 
streamer is simply connected and Riemann's mapping 
theorem applies; therefore the mapping exists. Since the 
boundary moves, the mapping is time dependent: we de- 
note it with z = ft{u>) where u> parametrizes the interior 
of the unit circle |a->| < 1. The complete map can be com- 
posed from a conformal map ( — h t (uo) that deforms the 
unit disc continuously, followed by the inversion z = 
Since h t {uj) is conformal on the unit disc, it is analytical 
and has a single zero which we choose to be at w = 0. 
Therefore /t(w) = l/ht(oj) has a single pole a and is 
otherwise analytical. Rather than a pole decomposition 
[20], we choose a Laurent expansion for ft(io) 



oo 

x + iy = z = / t (w) = ^2 a k( t ) U}k - 
fe=-i 



(5) 



This expansion allows us to identify the exact dynami- 
cal solutions (10) below. Taking a_i(i) as a real positive 
number makes the mapping unique, again according to 
Riemann's mapping theorem. 

(Hi) Now the potential on the unit disc can be cal- 
culated explicitly. Since /t(w) is a conformal mapping, 
the function &(z) is analytical if and only if the function 
<&(w) = $(/ t (w)) is analytical. The asymptote of $(lj) 
for u — > is determined by (2) and (5): for \x\, \z\ — > oo, 
we have ip(x,y) — ► —Eqx, hence $(z) — > —Eqz, and 
therefore with (5): $(lu) — > —E a-i(t)/cj for u — > 0. 
This means that the pole of $(w) at the origin of the 
unit disc u> = corresponds to the dipole of $(z) at 
z —* ±00. This dipole generates the field and the inter- 
facial motion. In the remainder of the unit disc, there 
are no sources or sinks of potential, hence $ is analytical 
there. Furthermore, at the boundary of the streamer, we 
have ip = from (3) or Re $ = 0, resp. The boundary of 
the streamer maps onto the unit circle, so Rc = 

for \u>\ = 1. Using the asymptotics at w —* and an- 
alyticity in the remaining region, the unique and exact 
solution for the potential is 



$(w) = E a-i(t) ( ui 



(6) 



(iv) The velocity Vbound. = cV<p (4) determines the mo- 
tion of the interface. This interface is the time dependent 
map ft (oj) of the unit circle u — e la parametrized by the 
angle a € [0, 27r). Therefore Eq. (4) determines the dy- 
namics of ft (e M ). According to [18,20], it is 



Rc 



id a f t *(e ia ) d t f t (e ia ) 



= c Rc 



(7) 
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The problem (l)-(4) is symmetric under reflection on 
the x-axis. We assume the solutions to have the same 
symmetry. Therefore all a k (t) have to be real. The po- 
sition (x,y)(a,t) of the point of the interface labelled 
by the angle a at time t can be read directly from the 
Laurent expansion (5) by inserting uj — e la ; it has essen- 
tially the form of a Fourier expansion of the unit circle 
where the circle and its position are created by the modes 
k = — 1 and 0: 

oc oo 

x(a, t) = ^ ctfe(i) cosfca , y(a,t) = ^ dfe(i) sinfca, 

k=-l fc=-l 

a k (t) real , a_i(i) > , < a < 2ir . (8) 

Substituting the mapping function (5) and the poten- 
tial (6) into the equation of motion for the mapping (7), 
and assuming the a k (t) to be real, we obtain for the evo- 
lution of the amplitudes a k (t) : 

oo 

^ k'a k '(t) d t a k (t) cos ((fc — k')a) 
fc,fc'=-i 

= 2£ co_i(t) cosa. (9) 

A closer investigation shows that this equation has 
an important property: suppose that the streamer 
boundary can be written initially as a finite series 

EfeL-i ffl fc(°) e4fcQ > a ^(o) ^ °- Then at a11 times t -> the 

interface is described by the same finite number of modes 

N 

z{a,t) = ^ a k (t) e ika , (10) 
fe=-i 

i.e., the ak(t) with k > N stay identical to zero at all 
times t > 0. Sorting the terms in (9) by coefficients of 
cos ka, the equation can be recast into N + 2 ordinary 
differential equations for the N + 2 functions a k (t) 

N-m 

^2 [i k + m ) a k+rn d t a k + ka k d t a k+m 
fc=-i 

= 2E c a_i S mA form = 0,...,iV+l, (11) 

where <5 TOi i is the Kronecker symbol. Eq. (11) is equiv- 
alent to a matrix equation of the form A({a k (t)}) ■ 

ft(a_i(t),...,Oiv(t)) = (0,2£ o ca-i(t),0,...,0), where 
the matrix A depends linearly on the {a k (t)}. 

Eqs. (10) and (11) identify large classes of analytical 
solutions with arbitrary fixed N. These solutions reduce 
the dynamical moving boundary problem in two spatial 
dimensions of Eqs. (l)-(4) exactly to a finite set of ordi- 
nary differential equations for the nonlinear coupling of 
the amplitudes a k (t) of modes e lka , < a < 2n. These 
equations are easy to integrate numerically or for small 
N even analytically. We will use this form to discuss now 
generic solutions of Eqs. (l)-(4) as the simplest approx- 
imation of a streamer in the Lozansky-Firsov limit. 



First, it is now easy to reproduce the uniformly pro- 
pagating ellipse solutions of [23,27] as the solutions with 
N = 1: for |ai| ^ |a_i|, the equations reduce to 9 t a_i = 
= dtCL\ and dta$ — 2Eqc a_i/(ai — a_i). These so- 
lutions correspond to ellipses whose principal radii are 
oriented along the axes. These radii maintain their val- 
ues r Xt y = a_i ± ai (assuming a_i > a\ > 0) and move 
with constant velocity i> e iiipse = —Eqc (r x + r y )/r y . The 
Lozansky-Firsov-parabola can be understood as limit 
cases of such uniformly propagating ellipses. 

In contrast to N < 1, all solutions with N > 2 have 
nontrivial dynamics. It can be tracked by integrating the 
N+2 ordinary differential equations (11) numerically and 
then plotting the boundaries (10) at consecutive times. 
Examples of such dynamics are shown in the figures. 
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FIG. 1. Upper panel: evolution of the interface in equal 
time steps up to time t = 0.1/ (Eoc) with initial condition 

a) z (a, 0) = e- ia + 0.6 ■ e M - 0.08 ■ e 2ia , 

b) z(a, 0) = z (a, 0) - 5 • 10~ 3 • e 8lQ 

c) z(a, 0) = z (a, 0) + 3 ■ 10~ 3 • e Sia 

d) z(a,0) = z (a,0) -4.5 • 10~ 7 • e 30ia , 

and lower panel: zoom into the unstable head of Fig. d. 

Fig. 1 shows four cases of the upward motion of a con- 
ically shaped streamer in equal time steps. The initial 
conditions are almost identical. On the leftmost figure, 
an ellipse is corrected only by a mode e 2la to create the 
conical shape. This shape with N — 2 eventually de- 
velops a concave tip, but only after much longer times 
than shown in the figure. In the other figures this coni- 
cal shape is perturbed initially by a minor perturbation 
with wavenumber 8 or 30, corresponding to N = 8 and 
30 in (10) and (11). The amplitude of the perturbation 
is chosen such that a cusp develops at time 0.1/(£bc). 
Depending on the sign of the amplitude, the cusp devel- 
ops on or off axis — where we stress that we are not 
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interested in the cusp itself, but in the earlier stages of 
evolution. Note that our reduction of the moving bound- 
ary problem to the set of ordinary differential equations 
(11) assures that the evolving shape is a true solution of 
the problem (l)-(4). Figs, lb and Id demonstrate that 
spontaneous branching is a possible solution. 
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FIG. 2. Evolution of the tip of an elongated "streamer" in 
equal time steps up to time t = 0.1/(Eoc); initial condition 
z(a, 0) = e- la + 0.9 • e la - 0.03 • e 2ia - 1.2 • 10" s ■ e 12la . 

In Fig. 2 the ionized body is longer stretched and only 
the tip is shown, again at 6 equidistant time steps. The 
streamer becomes slower when the head becomes flatter, 
since the electric field then diminishes. Eventually, the 
head becomes concave and "branches" . 

In summary, the solutions of the moving boundary 
problem (l)-(4) demonstrate the onset of branching 
within a purely deterministic model. They show a high 
sensitivity to minor deviations of the initial conditions. 
A streamer in the Lozansky-Firsov-limit is therefore also 
very sensitive to physical perturbations during the evo- 
lution, and simulations in this limit are just as sensitive 
to small numerical errors. But perturbations during the 
evolution are not necessary for branching. 

Our analysis applies to streamers in the Lozansky- 
Firsov-limit, i.e., to almost equipotential streamers that 
are surrounded by a very thin electrical screening layer. 
This limit is approached in our previous simulations [5,6]. 

These results raise the following questions that are 
presently under investigation: 1) When does a streamer 
reach this Lozansky-Firsov-limit that then gencrically 
leads to branching? 2) The formation of cusps should 
be suppressed by some microscopic stabilization mecha- 
nism. Is the electric screening length discussed in [5,30] 
sufficient to supply this mechanism? 3) If this stabiliza- 
tion is taken into account, can an interfacial model repro- 
duce numerical and physical streamer branching quanti- 
tavely? 4) How can the motion of the back end of the 
streamer be modelled appropriately (rather than assum- 
ing the velocity law v oc V<p (4) everywhere)? How can 
it be incorporated into the present analysis? 
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